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ANNA KOROLKO 

Abstract. We consider H (eisenberg)-type groups whose law of left translation gives rise to 
a bracket generating distribution of step 2. In the contrast with sub-Riemannian studies we 
furnish the horizontal distribution with a nondegenerate indefinite metric of arbitrary index 
and investigate the problem concerning causal geodesies on underlying manifolds. The exact 
formulae for geodesies are obtained. 



1. Introduction 

Sub-semi-Riemannian geometry is a study of nonholonomic manifolds M whose non-integrable 
tangent subbundles V of the tangent bundle TM are equipped with a nondegenerate indefinite 
metric g (sub-semi- Riemannian metric). If the metric g on T> is a positively definite metric, 
then the triple (M,T>,g) is known as sub-Riemannian manifold which has become an actively 
studied mathematical topic. If g is a nondegenerate indefinite metric on T> of index p G N, 
then the corresponding object can be called a sub-semi- Riemannian manifold; if p = 1, the 
manifold is called sub-Lorentzian. If D coincides with the entire tangent bundle TM and g 
is nondegenerate metric of index greater than 1, then the corresponding triplet reduces to 
semi-Riemannian manifold. There are numerous applications of nonholonomic manifolds with 
positively definite metric in problems in physics of high energy, control theory, electromag- 
netism, mechanical and biological systems with nonholonomic constraints. Therefore, it is 
interesting to start studying more generalized objects - sub-semi-Riemannian manifolds. 

Most of notions in sub-semi-Riemannian geometry originate in sub-Riemannian geometry 
\17\ [3] and Lorentzian geometry [18] . 

The intention of the article is to enlighten some aspects of sub-semi-Riemannian geom- 
etry. We start with a study of geometry, analysis and underlying algebraic structures of 
i/(eisenberg)-type groups with nondegenerate metric. These 2-step nilpotent Lie groups are 
interesting and natural examples of nonholonomic manifolds in sub-semi-Riemannian sense and 
at the same time first large class of objects with similar behavior which can be described in 
sub-semi-Riemannian geometry. The notion of H-type groups with positively definite metric 
was first introduced by A. Kaplan in his works [HI \TU[ [TT] with the help of theory of Clif- 
ford algebras. Such groups with a Riemannian metric provide nice examples of nonholonomic 
manifolds to study for a sub-Riemannian geometry [H [5] . 

Some examples of ii-type groups with semi-Riemannian metric on the distribution were 
already considered in [8j \12\ \T3[ [2] and this is the first attempt to obtain a generalizing result 
concerning such geometric objects of interest as geodesies. We solve the geodesic equation for 
sub-semi-Riemannian geodesies, which coincides with the geodesic equation in sub-Riemannian 
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case if the considered metric tensor is Riemannian. The essential role here plays the skew- 
symmetry of emerging objects with respect to the semi-Riemannian product on V, which 
is not positively definite anymore (see the section [4]) . The skew-symmetry arises from the 
use of Clifford algebras and their representations. A good deal of the article is devoted to 
the description of skew-symmetric transformations w. r. t. semi-Riemannian metric. The 
canonical view and, therefore, eigenvalues of corresponding matrices that are skew-symmetric 
w. r. t. semi-Riemannian metric g greatly depends on the index of g. 

The article is organized in the following way. In Section [2] we introduce H-type groups 
endowed with semi-Riemannian metric. We also present necessary basic notions from space- 
time geometry. In Section [3] we indroduce exponential map, with its help we calculate left- 
invariant vector fields and the connection. Section [4] is devoted to calculations of geodesies and 
auxiliary algebraic backgrounds. And in Section [5] we provide the example of -ff-type group 
with semi-Riemannian metric which illustrates how the expounded theory could be applied for 
the problem of finding geodesies in a particular case. We conclude the survey by the list of 
open problems in Section 6. 



2. General structure of Heisenberg-type groups with nondegenerate metric 

Let U and V be vector spaces of dimensions m and n respectively with positively definite 
product {-,-)u on U and nondegenerate product {-,-)v of signature (p, q), p + q = n,onV. 
Denote by r\ the metric tensor of the product on V, which is a diagonal matrix with first p 
negative unities and then q positive unities, i. e. 



(2.1) r, 



-I P 
I q 



where I p and I q are usual identity matrices of dimensions (p x p) and (q x q) respectively. The 
corresponding metric g : V x V — > M is called a nondegenerate metric of index p, p £ N. Then 
as usually g(v,w) = (v,w)v = w t r]v, where v,w 6 V and w l denotes the vector transposed to 
w. Nondegeneracy is understood here in the sense that if {v,w)y = for any w 6 V, then 
v = 0. If the index of the metric is 1, then it is called the Lorentzian metric. The metric with 
the index is just a Riemannian metric. In the case of arbitrary index such a nondegenerate 
indefinite metric received the name semi-Riemannian metric. 

Let us denote by the symbol End(V) the set of linear maps from V to V. Consider a linear 
mapping j : U —> End(V), which satisfies the following properties: 

^ 22 ^ 3 2 {u) = -\u\ulv, 

ii) (i]j(u)v,w) v = -(v,i]j(u)w) v , 

where u £ U, v, w £ V and Jy is a (n x n) identity matrix in V, \u\fj = (u,u)u- Conditions 



(2.2) are taken from |11| and adapted for the case of semi-Riemannian metric on V. For all 
u G U the linear map j(u) provides a complex structure on V compatible with its euclidean 
structure. By abuse of notations we use the symbol j to denote both the mapping and the 
corresponding (n x n)-matrix. 

Remark. The property of j being skew-symmetric remains the same as in the case of H- 
type groups with positively definite metrics on both spaces U and V ([H]). Indeed, notice 
that 

(r]j(u)v,w)v = w t r]rjj(u)v = w t j(u)v, 
(v,r)j(u)w)y = w t j(u) t r] t r]v = w t j(u) t v 
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and from the condition ii) it follows that 

(2.3) j(u) = -j( u y. 

Then we can define the bilinar skew-symmetric map [-,•]: V x V — > U by 

(2.4) (u, [v, w])u = {rjj(u)v, w)y. 
We can now define a Lie algebra 

Q = U®V 

with a Lie structure on it: 

[u + v, r + w] = [v, w] (u, r G U, v, w G V). 

Constructed in this way Lie algebra is graded and 2-step nilpotent with the center U. As in 
we shall refer to these algebras as -ff-type algebras with nondegenerate metric (•, -)y. 
Recall that if U is a m-dimensional real space M. m with a basis consisting of m orthogonal unit 
vectors e%, . . . , e m , then the Clifford algebra C1(C7, m) is the real associative algebra generated 
by the set {ei, . . . , e m }, satisfying the relations (see [16J) 



e 2 a = — 1, a = 1, . . . , m, 

Ca^p = — ege a , for a ^ (3. 

The Clifford algebra Cl(U, m) is 2 m -dimensional with the following basis: 

1 the scalar 
ex,...,e m vectors 

eie 2 , • • • , e m _ie m bivectors 

eie 2 e 3 , . . . , e m _ 2 e m _ie m 3-vectors 

eie 2 . . . e m the m- volume element. 

Hence an arbitrary element of the basis may be written as e_4 = e ai . . . e ai here A = 
(«i, . . . ,Oj) and 1 ^ a.\ < a 2 < . . . < oti ^ m. All elements of this type share the basic 
property of the imaginary unit of being the square root of —1, i. e. (e^) 2 = —1. It is known 
that real Clifford algebras C1(C/, m) are isomorphic to real matrix algebras with entries in R, 
C, H or in IR © M, H © H, where by K, C and EI we denote division rings of real numbers, 
complex numbers and of quaternions respectively. 



The mapping j : U — > End(y) satisfying (2.2 ) extends to a representation of Clifford algebra 



Cl(U, m) on V; in other words, the space V becomes a Clifford module over Cl(U, m). The 
images of the generators e±,. . . ,e m of a Clifford algebra Cl(U,m) under a representation j: 
ja '■= j(e a ), a = l,...,m, are called (standard) generators of the representation. It was 
noticed that without loss of generality one can represent standard generators as a set of anti- 
commuting skew-symmetric matrices with squares —ly (see [2]). Moreover, these matrices can 
be expressed as homogeneous tensor products of the standard Pauli (2 x 2)-matrices, though 
the determination of such matrices is nontrivial. 

Recall also that V admits a structure of Clifford module over C\(U, m) if and only if the 
following condition on dimensions n and m of vector spaces V and U, respectively, holds: 
^ m < p(n), where p is the Hurwitz-Radon function defined by n = k2 4r+s i— > p(n) = 8r + 2 s , 
k - odd, ^ s ^ 3 [11]. If V is a normed division algebra, then U = ImV is a space of imaginary 
elements of V and dim U = m = n — 1. Ifm = and n is arbitrary, then Q is the Euclidean 
n-dimensional space; if m = 1 and n is even, then Q is the Heisenberg algebra; if m = 3 and 
n = 4k, then Q is the Quaternion algebra; if m = 7 and n = 8, then we call the corresponding 
algebra an H © H-algebra. These four algebras generate a special class among the .ff-type 
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algebras, satisfying the j 2 -condition: for all ui,u 2 G U with (ui, 112)1; = and all nonzero 
oeF there exists a U3 G U so that 

i(^l)j'( U 2)^ =j(u 3 )v. 

Exact formulae for geodesies on some of corresponding groups with semi-Riemannian metric 
on V were found in \12\ [TBI I14| . In the present article we concentrate our attention on the 
general .H-type groups with nondegenerate metric on V. 

Let us consider the H-type group G with the Lie algebra Q = V © U and nondegenerate 
product (•, -)v on V . The Lie algebra Q is identified with the tangent space to G at the identity 
T e G. The scalar product (•, •) on the algebra Q is a sum of the inner product (•, -)u on U and 
nondegenerate product {-,-)v on V. Push-forward allows to define the scalar product on the 
whole tangent bundle TG = yj a T a G, not only in T e G: 

(;-)(a) = (dL a -i;dL a -i-)(e), 

where L a denotes the left translation on G by the element a G G. Since V C T e G, the left 
translation of V by element a is dL a (V) C T a G. The mapping V: a h-> dL CT (V) is called a 
horizontal distribution and any vector -y(cr) G dL a (V) is called a horizontal vector at point o\ 

By analogy with sub-Riemannian geometry the absolutely continuous curve c(s) : [0, 1] — > T> 
is called the horizontal curve if the tangent vector c(s) belongs to T> at each s G [0, 1]. The 
Lie group G was introduced in such a way, that the horizontal distribution T> possesses a 
2-step bracket generating property. It means that for any point p G G we have that any 
tangent vector can be presented as a linear combination of the vectors of the following types 
v(p), [v,w](p) G T p M, where all vectors v, w are horizontal. A result of Chow [6] says that 
any two points on G can be connected by a piecewise smooth horizontal curve because of 
the bracket generating property of T>. Chow theorem guarantees that the connectivity of a 
manifold by curves tangent to a given distribution T> does not depend neither on metric defined 
on it nor on tangent bundle. 

In this article we study the triple (G, T>, g), which we call the sub-semi-Riemannian manifold, 
where the nondegenerate bundle-type metric g of arbitrary index p is given on distribution T>. 
We will call the corresponding group G sub-semi Riemannian H-type group. Due to left 
translation we can work further with the distribution at the point e, i. e. with dL e (V) = V 
instead of T>. 

The nondegeneracy of the metric g brings the causal character of the vectors in V . The 
notations are carried from the Lorentzian geometry [dj. Even though the metric g has the 
arbitrary index p and the time orientation can not be defined in general (unless p = 1), we 
can anyway introduce the notions of timelikeness, spacelikeness and lightlikeness. A horizontal 
vector v G V is called timelike if (v,v)y < 0, spacelike if (v,v)y > or v = 0, lightlike if 
(v,v)v = and u^O, nonspacelike if (v,v)y ^ 0. A horizontal curve is called timelike if its 
tangent vector is timelike at each point; similarly for spacelike, lightlike and nonspacelike. 

3. Exponential map, left-invariant vector fields and connection 

We define the exponential map as a function exp : Q — > G in a usual way in the Lie groups 
theory. Then we can define smooth maps u : G — > U and v : G — > V such that 

a = exp(u(<r) + v(<r)) for a G G. 
Since Q is a 2-step nilpotent Lie algebra, then the Baker-Campbell-Hausdorff formula yields 
v(cr 1 cr 2 ) = v(cji) + v(er 2 ), 
(3.1) 1 

u(o-i<72) = 11(01) + U(0"2) + 7T[v(0l), v(<T 2 )]. 
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Suppose that {d Vj }™ =1 and {d Ua }™ =1 are orthonormal bases for V and U respectively such 

that (d Vi ,d Vj ) v = £jSij and (d Ua ,d Ufj )u = 8 a p, (d Vi ,d Ua ) = 0, <% is a usual Kronecker symbol 
and £j = — 1 for 1 ^ j ; ^ p and £j = 1 for p+l^j^p + q = n. Denote by Afj the analog of 
Clifford coefficients, i. e. the coefficients in the decomposition 

n 

Denote by the structure constants in 

m 

'd Vi ,d Vj ]=j2K 9 ^ eU - 

(3=1 

They remain the same as in case of H-type groups with sub-Riemannian metric. The structure 
constants B^- are connected to the analog of Clifford coefficients in a following way: 

Lemma 1. 

R a — p • A a 

ij ~ fc J yl «j' 

Proof. We calculate 

m ra 

(3.2) {m(d Ua )d Vi ,d Vj ) v = (d Ua , [d Vi ,d Vj ])u = Y,B^(d Ua ,d Ufj )u = J>g<W = B ?v 

(3=1 p=l 

but from the other side 

n n n 

(vj{du a )d Vi ,d Vj ) v = {J2 A ?kdv k ,d Vj ) v = ^A%(d Vkl d V] ) v = Y^A%e k 5 kj = e 3 A% 



k=l k=l k=l 



Let us denote the global coordinates on the group G by {v±, . . . , v n , u±, . . . , u m }: 

n m 

v(o) = ^Vid Vi {e), u(cj) = ^u a d Ua {e). 



i=l a=l 

Let us compute the Lie bracket 



n n 

[v(ai),v(<7 2 )] = \^v}d Vi {e)^v 2 j d Vj {e) 

i=l j=l 

n n n 

= ^ 2 ( E v ) B h d ^ + E + • • • + E v ) B T^ m ) 

j=i j=i j=i 

n n n 

+ v l ( E v ) B h^ + E + • • • + E v]B%d Um ) 



j=i j=i j=i 



n n n 

+«» ( E V ) B )n^ + E V ) B lndu 2 +-.. + E ^>A m ) ' 



□ 
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Therefore, the left-invariant vector fields obtained by the rule Vi (a) = dL a d Vi (e) are as follows 

- m n 

^i = ^+2E(E v 3 B< ji) d u a i i = l,...,n 

a=l j=l 

and 

U a = d Ua , a = l,...m. 

The vector fields remain the same as in sub-Riemannian case when the metric on V is positively 
definite, since the matrix dL a is obtained from the group structure, which we do not change 
(compare with [H]). Matrices B a = (Bfy, a = l,...,m, i,j = l,...,n, are in one-to-one 
correspondance defined by (3.2) with matrices j a = j(d Ua ) that are the standard generators 
of representation for Clifford algebra Cl(U,m). Analogously to sub-Riemannian case we can 
write the horizontal gradient, i. e. the vector V = (V\, . . . , V n ), in the form 

V = S7 v + ^B 1 vd Ul + ... + B m vd Uri 

with v = (vi,...,v n ) and V„ = (d Vl , . . . ,d Vn ). 

Remark. Notice that {V^, U a }, i = 1, . . . , n, a = 1, . . . , m is an orthonormal basis with 
respect to the left-invariant metric (•, •). 

Proof. By definition of metric (-, •) on algebra Q we have: 

(Vi,Vj) = (Vi,Vj) v = (d n ,d Vj ) v = EjSij, 
(U a ,Up) = (U a ,Ufs)u = (d Ua ,d U p)u = 5 a p, 
(Vi,U a ) = 0. 

□ 

The connection on a group G is defined by Koszul formula Q18J) 

(V*-y, Z) = l - (X{Y, Z) + Y(Z, X) - Z(X, Y) - (X, [Y, Z\) + (Y, [Z, X}} + (Z, [X, Y])). 
For the orthonormal basis of Q this formula receives the form 

( Vx Y Z) = ±(-{x, [Y z\) + (Y, [z, x]) + (z, [x, y]>) , 

where X, Y, Z are arbitrary basis elements in the corresponding Lie algebra Q = U ® V. In the 
case of H -type group we calculate that 

(WiVj,U a ) = \(U a ,[Vi,VA)v = ^(vj(U a )Vi,Yj) v , 
(Vu a V i ,V j )=-^(U a ,[V i ,V^) u = -^( v j{U a )V i ,V j ) v , 

(WiU^Vj) = \(U a , [Vj,Vi\)u = ~(U a , [Vi,Vj])u = -~(vj[U a )Vi,Vi)v, 

for any Vi,Vj 6 V, i / j, and U a G U. It is easy to verify that all the other combinations 
result as 0. From here we obtain that the connection is given by 

w«y;- = ~[^], 

(3 ' 3) Vu a ^ = S7v.Ua = ~qj[U a )Vi, 

Vu a Up = 0. 
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4. Geodesics 

In this section we compute geodesics passing through the identity with a given initial vector. 
By geodesic we mean a curve 7: [0, 1] — >• G which satisfies the geodesic equation V-y7 = 0- We 
define functions t 1— > v(t) and t *— > u(t), such that j(t) = (v(t),u(t)) and v(t) and u(t) are the 
projections of 7 (i) onto V and U respectively. 

Then in terms of global coordinates on G one gets j(t) = (t>i(i), . . . , v n (t),u\(t), . . . , u m (tj) 
and the derivative with respect to t is 



IV I I V IV ^ I I V IV I I V 

^ i=l a=l i=l a=l j=l a=l 



n m 1 n 



= iid vi + J2{ iia + \l2 r j''j'' r ) = v+(u+ - [v,v]) . 

i=l a=l *,i=l 

Therefore, 

(4.2) Vi7 = VjV + \7j(u + -[v,v]). 

Let us compute the first term in right hand side of (|4.2l) 



V-y* = S7i + u+i[v, v f = S7vV + X7uV + - V[*,«] v = v- ^j{u)v - ~i]j({v,v])v 

.. 1 .,. l r . n . 

= v - -VJ{u+ -[v,v\)v. 

Analogously, we calculate the second term 



Therefore, 



Vj(u + ^[^D = il +^[v,v] - -r)j(u + £ 



V-y7= (p-rij[u + -[v,v])vj + (u + ~[v,v}) 



and the geodesic 7 is characterised by the system of equations: 

v - r]j(u + hv,v])v = 0, 
(4.3) 1 2 

u + -[v,v] = 0. 

From the second equation of this sytem we can deduce that u + \ [v, v] =constant. It is enough 
to describe geodesics which pass through 0, since the left translation of every geodesic is again 
geodesic. Let us assume the following initial conditions for 7: 

7(0) = 0, 7(0) = v° + u°, 



which implies that 



Then 



u(0) = 0, «(0) = 0, 
v(0) = v°, u(Q) = u°. 



(4.4) ii+ -[v,v] = n° for all t G [0, 1] 
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and the first equation of (4.3) takes the form 
(4.5) v — 7]j(u )v 



0. 



Observe that this equation coincides with the one for sub-Riemannian geodesies when r\ is an 
identity matrix in V 



To find a solution v(t) of (4.5) let us first recall some facts from linear algebra. 



4.1. Skew-symmetry with respect to semi-Riemannian metric. Recall that j(u) can 
be written as a real skew-symmetric (n x n) matrix: j = —y if we fix a basis on V. Therefore, 
it has purely imaginary eigenvalues ±ifJ>ki Mfe € R, k = 1, . . . , ^, since n is an even number. It 
is possible to bring every skew-symmetric matrix to a block diagonal form by an orthogonal 
transformation. Thus, j can be written in the form 



(4.6) 



J 



QjQ\ 



where j is a Murnaghan's canonical form of j: 

I mi 
-Ml o 



J 







V 











112 
-H2 










- fin/2 



/ 



and Q is an orthogonal transformation. From the condition j 
3 



-J 



we obtain also that 



-j . And since j(u) has the property i) of (2.2), all the eigenvalues of j satisfy the 



condition \i k 



! ,k = l,... 



n 
' 2 • 



Denote by A the matrix rjj(u) from (4.5). The condition ii) of (2.2) yields that the matrix 



A is skew-symmetric with respect to the scalar product on V and it can be rewritten as 
w l r]Av + (AwYrjv = 0, W, w £ V, from where we obtain that 

(4.7) rjA + A l r] = 0, 

that is, A t = —r\Ar\ (skew-symmetric w. r. t. rj). 

Every square matrix can be represented as a sum of symmetric and skew-symmetric ones, 



but in our case these matrices are of very special form, that follows from (2.3): 



(4.8) 



J 



p 











+ 











Ai +A 2 , 



where <P p and <P q are matrices of dimensions p x p and q x q respectively such that <P P 
and <P q = —<P q ; & is a (p x g)-matrix and, therefore, matrix A 2 is skew-symmetric. In analogous 



p 



way using (|4.7|) we may represent the matrix A 
(4.9) 



-I P 
o /„ 



p 

-xpt 



®P 



rjj as follows 











+ 





-1^* 





il +^2 



The essential distinction between matrices j and A = rjj is that in the decomposition (4.8) of 



j both matrices A\ and A 2 are skew-symmetric, but in decomposition (4.9) of A the matrix 
A\ is skew-symmetric, but A 2 is symmetric. 
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The condition (4.7) is the condition for A to be skew-symmetric matrix with respect to 



(v)y and is analogue of skew-symmetry condition (2.3) for positively definite metric. Using 



it now together with the condition i) of (2.2) we get 

j 2 {u) = rjAnA = -A l r\r\A = -A 1 A = -\u\ 2 I V - 

Prom another hand, 

f[u) = r]Ar)A = -r]AA t T] = -\u\ 2 Iv, 
from where it follows that AA l = r]\u\ 2 Ivn = |u| 2 /v/- Therefore, 
(4.10) AA % = A 1 A = \u\ 2 I v . 



This allows us to rewrite Kaplan's conditions ( 2.2 ) as follows. Consider a mapping A u : V —> 
V (we will omit index u for the sake of simplicity) , which is assigned to each u G U and satisfies 
the following properties: 

i') AA l = A 1 A = \u\ 2 I v , 
ii') (Av, w)v = — (v, Aw)v- 

The calculations above show that the conditions i) and i') are equivalent; the condition ii') is 
just a rewritten form of ii) in terms of matrix A. Therefore, we get an independent definition 
of sub-semi- Riemannian H- type group. 



Another effect of the property (4.10) is that A appears to have n mutually orthogonal unit 
eigenvectors and is a diagonalisable matrix (over C in general). Observe that the matrix A 2 is 
not symmetric, but it satisfies the condition {A 2 ) 1 = r]A 2 r] (symmetric with respect to rj). 

It is also worth mentioning that even though the matrix A is real, the condition A 1 = A* 
for adjoint matrix A* is not true in our situation as in the case of positively definite scalar 
product. We define the adjoint matrix A* ,r) in our case by (Av,w)y = (v, A*' TI 'w}v- It has the 
following view A*' ri = 7]A t 7j, i. e. A*' v = —A. Then the condition of A being normal holds: 
AA*' 7 ! = A*' ri A = -A 2 . 

Let us write the characteristic polynomial for the matrix A: 

(4.11) PA (X) = X n + fon-iA"" 1 + . . . + b 2 X 2 + h\ + b = 0, 

where the coefficients in front of X k , k = 0, . . . , n — 1, are computed by the known formula: 



n 

=i 



where by symbol A l k we denote (k x A;)-dimensional principal minor of the matrix A, that 
is the determinant of the matrix obtained from the matrix A by erasing its n — k diagonal 
elements. For a (n x n)-square matrix A there are n (k x k)- dimensional principal minors for 
each k = 1, . . . , n. Notice that if k is odd, n — k is also odd. 

Lemma 2. All odd- dimensional principal minors of skew-symmetric w. r. t. n matrix A are 
zero. 

Proof. We calculate 

det(4)* = det 4 = det ( - ^(4)^) = det ( - (4)*) = det(-I fc ) det(4)* = {-if det(4)*, 
where It is a (k x A;)-identity matrix. From here it follows that det A l k = for odd k. □ 
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Therefore, the characteristic polynomial (4.11) receives the following form: 



0, 



(4.12) PA (X) = X n + a f _xA n - z + . . . + aiA z + a 

where = b2k, k = 1, . . . , ^ — 1, and ao = 60 = det A. Observe that ao ^ since the dimension 
n of the matrix A is even. Therefore, is not an eigenvalue of A and all eigenvalues of A come 
in pairs ±Afc, k = 1, . . . \. 

Lemma 3. // A is an eigenvalue for A with eigenvector v, then —A is an eigenvalue for A 1 
with eigenvector r\v. 

Proof. Notice that (v,w)v = (vViV^V for any two vectors v and w in V. We calculate 

(Av, w)y = (Xv, w)y = X(v, w)v = X(rjv, rjw)v, Vu> 6 V. 
But from another hand, 

(Av,w)v = (—r]Ar] t v,w)v = — (A t rjv,r}w)y, \/w € V, 
yielding the desired result. □ 

The operator A has n in general complex eigenvalues and the corresponding complex eigen- 
vectors. The following fact is known, see, for example, |19j . 

Proposition 1. Let A be a (n x n) -matrix with real entries. If the equation det (A — XI) = 

has a complex solution A = a + i/3, f3 7^ 0, then there exist nonzero vectors x and y in V such 
that 

Ax = ax — (3y, Ay = ay + fix. 
In particular, spanjx, y} is a nontrivial invariant subspace for A. 

Notice that if A = a + i(3 is a solution of characteristic polynomial (4.12) with multiplicity 
I, then its conjugate A = a — i/3 is also solution of (4.12) with multiplicity Due to this fact, 



lemma [3j proposition [T] and the fact that A is diagonalisable over the complex space V, the 
following decomposition of A is possible: A = PDP~ l , where P is orthogonal matrix and D 
is a block-diagonal matrix with (2 x 2)-blocks of the following possible forms: 



(4.13) 
(4.14) 

and (4 x 4)-blocks of the form 
/ 

(4.15) 



A, 






-Xi 



-IV A 



l,...s, 



j = s + 1, . . .r, 



-fik a k 



\ 







-ak 

Pk 



-Pk 
-ak 



n 
4 



r 
2' 



where the eigenvalues are rearranged in the following order: 

• ±Ai, . . . , ±A S , for some s S {0, . . . , |}, are real eigenvalues of A, 



±A s+ i — ±iv s+ \, . . . , ±A r 
eigenvalues of A, and 



±iv r , for some r G {0, . . . , S}, r ^ s, are purely imaginary 



±A, 



r+l 



±(ai + ifii), ±A r+ i, . . ., ±Ai 



), ±A«_i: are complex 



eigenvalues of A corresponding to 



blocks of the dimension (4x4). 
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Notice that the integer r is necessarily even. It is natural that (4 x 4)-dimensional blocks 
associated with complex eigenvalues ±(a±i/3) can appear only if the dimension n is bigger or 
equal to 4. 

Observe that for real eigenvalues Aj the matrix ( ' \ J , i = 1, . . . , s, is a diagonalised 



form for the matrix 





A; 



A, 




Wj 





-IV A 



is a diagonalised form of 



x ~ Xi , 

; and for purely imaginary eigenvalues Xj 

v, 



ivj the matrix 







j = s + l,...,r. 



Notice that even though A is skew-symmetric w. r. t. 77, i. e. the condition (4.7) holds, the 



matrix D does not have this property in general as long as 77 is not identity matrix. 
Lemma 4. All eigenvalues of A have absolute value \u\. 

Proof. Let us calculate the product AA t = P D P~ l (P l )~ l D l P l = PDD t P t . Observe that the 
matrix DD t is a diagonal matrix with entries |A* 



k = 1, . . . , \ . Then, using the property 
(4.10|), we see that all eigenvalues of A, real and complex ones, have the same absolute values: 
u| 2 , k = !,...,§. □ 



Taking into account previous comments about eigenvalues we can write that 

\ 1 = . . . = A s = v s+ i = . . . = v r = \u\, 
s,r € {0, . . . , § } and a\ + /3? = |u| 2 , k 



1 - 



r 
2' 



Thus, we showed that the matrix D is built from the the following cells: 





u 



u 




u 




or 





V 



a /3 
-j3 a 








—a 

P 



\ 


-a J 



where a 2 + /3 2 = |-u| 2 . 

The appearance of one or another cell of these types in the matrix D depends on the index 
p of the metric g. 

Below we show that the number s of pairs of real eigenvalues of A is equal to 1 in the case of 
odd index p of the metric 77, and s = in case of p is even. In other words, if index p is even, 
then there are no real eigenvalues of A, and if p is odd, there are only two real eigenvalues that 
are of the form ±|n|. 



Note. Notice that det^4 = det D 



\u\r(\u\r~ s (m 



= (-l) s M n . But on the 

other hand from (4.6) we have det A = det (r]j) = det (rjQjQ -1 ) = (— l) p |n| n . Thus, if p is 



odd, then s is odd too, and if p is even, then s is also even. Notice also that r — s has the same 
parity as s. 

Remark. We restrict our consideration only to the case when the index of the metric p is 
less or equal to ^, because in an opposite situation we can study the case with the metric —77, 
which is analogous. 

4.2. Matrices ji, . . . , j m - We remind that matrix j(u) does not depend on the metric r\ and 
is the same as for sub-Riemannian //-type group, i. e. can be written via endomorphisms 
jli ■ ■ ■ ijm represented by structure constants ( |3.2[ ) in the following form: 

u = (u 1: . . .,u m ) € U. 



J( u ) 



U\3\ + . . . + U m J r 



12 ANNA KOROLKO 

Each of j a is, what is called in linear algebra, a one-column and one-row matrix with 1 at the 
intersection of i-th row and j-th column, —1 at the intersection of j'-th row and i-th column, for 
i 7^ j and each i,j < n and zeros everywhere else (see, for example, matrices j a , a = 1, . . . , 7 
from Section [5]); and satisfies the following properties 

(4.16) J* = -j a , jl = -Iy, j a jp = -jpj a , a^fi. 

The block-diagonal (n x n)-matrix consisting of blocks of type ( 1 n ] is one of them, 



and let us denote it by j\ for simplicity. Since Cl(U, m) by definition is the associative algebra 
freely generated by m anti-commuting square roots of minus unity, one can associate matrices 
jii ■ ■ ■ i3m with square roots of — Iy. Their linear combination represents the corresponding 
endomorphism j(u). The matrices j a have n eigenvalues which come in pairs ±i and, therefore, 
det j a = 1 for each a = 1, . . . , m. 

We remind that j(u) is a (n x n)-matrix with the following properties: 

• f(u) = -j(u); 

• = -\u\ 2 Iv, 

• the eigenvalues of j(u) are all of the form ±i|u|; 

• the determinant of j(u) det j(u) = \u\ n . 

From here we can catch out the following properties of matrices rjj a , a = 1, . . . , m. 

Lemma 5. Each matrix r\j a is skew- symmetric w. r. t. rj, i. e. rjj a + j a rj = 0. 



Proof. The proof is the same as for an analogous statement (4.7) for the matrix A = r}j{u) in 
the previous sub-section. □ 

Analogous argument as for matrix A can be applied to show that all eigenvalues of matrices 
VJa, a = 1, . . . , m, have absolute value 1 and come in pairs ±A^, k = 1, . . . , ^. 

Lemma 6. All eigenvalues for matrix r]j a are of the form ±1 or ±i, for each a = 1, . . . , m. 

Proof. If a + ib is an eigenvalue for r/j a , then by proposition [T] there exist vectors x and y such 
that rjj a x = ax — by and rjj a y = ay + bx. Thus, on the space spanned by vectors x and y the 
linear transformation r/j a has the matrix 

a b 
—b a 

which is possible only if a = or b = 0, since r]j a is a 1-row and 1-column matrix. Due to the 
fact that all eigenvalues of rjj a have absolute value 1, they should be of the form ±1 or ±i. □ 

Lemma 7. Matrices r]j a , a = 1, ... ,m, anti-commute with each other. 

Proof. Using lemma [5] and anti-commutativity of matrices j a we calculate for a 7^ (3: 



□ 



Lemma 8. // the index p of metric 77 is odd, then each matrix r]j a has at least two real 
eigenvalues ±1. 
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Proof. The proof follows from the fact that det = (— l) p = — 1 for the odd p. From 

another hand, the determinant of matrix nj a is equal to product of all eigenvalues of rjj a that 
is (— l 2 ) • . . . • (— l 2 ) (— i 2 ) • . . . • (— i 2 ), where ^ s ^ S is an integer. This product equals to 

s times ^—s times 

— 1 if and only if there is an odd number s of pairs of real eigenvalues dbl. □ 

Note that matrix rjjx has only 2 real eigenvalues ±1 if p is odd. For p - even there exists at 
least one matrix from rjji, . . . , rjj m which does not have real eigenvalues, for example, matrix 
VJi- 

Notice also that if the matrix nj a has eigenvalues of the form ±1 and/or ±i, then the matrix 
UaVJa has eigenvalues of the form in„ and/or ±iu a respectively, where u a is a real number. 

Lemma 9. 7/±l are eigenvalues for all of the matrices nji, . . . ,i]j m , then ±\u\ are eigenvalues 
for r,j(u). 

Proof. We will need the following proposition. 

Proposition 2. Given two anti- commuting (n x n)-matrices A and B with real eigenvalues 
±a and ±6 respectively, their sum A + B has two of its n eigenvalues of the form ±V« 2 + b 2 . 

Proof of the proposition. Let V a be the a-eigenspace of the matrix A. Then by anti- 
commutativity we find BV a Q V- a - Notice that matrices A and B 2 commute and, therefore, 
are simultaneously diagonalizable. So there exists an eigenvector v 6 V a with the property 
B 2 v = b 2 v. Then notice that for the vector w = tBv we have the following relations 

Av = av, Bv = bw, 

Aw = —aw, Bw = bv. 

From here it follows that on the space spanned by vectors v and w the linear transformation 
A + B has the matrix 

a b 
b —a 

which has eigenvectors with eigenvalues ±V a 2 + b 2 . □ 
Easy induction shows that the same argument applies to the sum of m matrices: if {^ajJJLi 

is a set of anti-commuting matrices and each matrix A a has two real eigenvalues ±a Q , a = 

1, . . . , m, then the matrix A\ + . . . + A m has two real eigenvalues of the form ±\/ a 2 + . . . + a^. 
We obtain the statement of Lemma [9] by taking matrices nj a = A a and numbers u a = a a , 

a = 1, . . . , m. □ 

Thus, Lemmas 5]|9 can be summed up in the following theorem. 

Theorem 1. If the index p of the metric n is odd, then there are two real eigenvalues ±|ii| for 
the matrix rjj(u). If p is even, then there are no real eigenvalues for rjj(u) at all. 

Proof. In order to have 2 real eigenvalues ±|u| for the matrix A = nj(u) all matrices rjj a should 
have ±1 as eigenvalues, which is possible only if p is odd. There are only two real eigenvalues 
for A, since matrix rjji has only two real eigenvalues. If p is even, matrix rjji does not have 
real eigenvalues, and, therefore, A does not have real eigenvalues either. □ 

Further idea is to construct 5 2-dimensional invariant subspaces of the operator A: V — > V . 
Remind that V has a basis consisting of p timelike and q spacelike vectors. There are two 
possible cases here depending on the index p being odd or even. In case p is odd we have 

v = w e Wi e . . . e wvi e w x e . . . e wvi , 

2 2 
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where the restriction of the metric g on the subspace Wo g\ Wn is Lorentzian, that is it has 
index 1, g 



p-l 
' 2 



, is a negatively definite metric and g\^r., 3 



1 , — , , is a 



positively definite metric. In case of even p we have 

V = Wi@ ...@Wz@Wi 



Wi. 

2 



where g\ w , i = 1, . . . is a, negatively definite metric and g\y^., j = 1, . . . , |, is a positively 
definite metric. The basis of the subspace W$ consists of one timelike and one spacelike unit 
vectors, the basis of subspaces Wi, i / 0, consists of timelike unit vectors and the basis of 
subspaces Wj consists of spacelike unit vectors. 



4.3. Formulae for geodesies. Let us return to the system of equations (4.3). It is equivalent 
to the system in new basis v = Pv £ M. n 



(4.17) 



Dv, 



where D = P~ l AP is a block-type matrix with blocks of the form (|4.13l)-(|4.15l) and P is 
orthogonal matrix. 



It is convinient to remunerate coordinates in v 



in the following way: 



l)n ir 1)5 if VTL-L.L ,VBl-1-L 
4 T 21 4 T 22 4 T 23 4 T 24 



V = [V S1 ,V S2 ,V s+ l 1 ,V s+ l 2 , . . . , V ri , V r2 , V r+ll , V r+ l 2 , V r+la , V T+U , 

where first coordinates v Sl ,v S2 correspond to real eigenvalues ±|u| (if any), then coordinates 
Vji > Vj 2 , j = s + 1, . . . , r, corresponding to purely imaginary eigenvalues ±i\u\ (if any), and then 
coordinates v r+kl , v r+k2 , v r+ks , v r+ki , k = 1, . . . , § 



| corresponding to complex eigenvalues 
+(a k ± i/3 k ), a\ + fi\ = |u| 2 , k = 1, . . . f - § (if | < |). Here lowest index in v ix , V{ 2 , 6j 3 , 
denotes the number of the coordinate in v associated with the eigenvalues Aj, — Aj, Aj, — Aj 
respectively. 



Therefore, the system (4.17) splits into S independent systems each in 2 variables 



(4.18) 



V S1 = 








Vs 2 = 


\u\v si , 






v h = 


\u\Vj 2 , 






Via = 


-\u\v h , 


3=8+1,.. 


■,r, 




= OikV r+kl 


+ PkVr+k 2 > 




«r+fc 2 


= —Pk^r+ki 


+ akV r+ k 2 , 






= —a k v r+k3 


— PkVr+k^, 


al + Pl = \u\ 2 , 




= fikV r +k 3 


— Ct k V r+ k 4 , 


h — 1 n r 
^ — L i ■ ■ ■ > 4 2 



There are n equations and n variables in (4.18 ). We are looking for the solution v = v(t) G V, 

5 . From the system (4.18) we get 



t G [— oo, +oo], v(0) = and with initial velocity v(0) 



v 



J S 2 



-c Sl e-H* + c e \u\t ~o 



-\u\t 



+ c S2 e 



\u\t 



"S 2 



C S\ + C S2 , 
C Sl + C 



Vj 1 = Cj sin(|-u|i) + Cj 2 cos(|n|t), 
vj 2 = cj cos(\u\t) — Cj 2 sin(|u|t), v 



*2 ' 



5° 

n 

- c 
j2 °n ' 

s + 1, . . . ,r, 
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e akt ( - c r+kl cos fat + c r+fc2 sin j3 k t) , 
e akt (c r+kl sin f3 k t + c r+k2 cos /3 fc t) , 



It follows that 



e Qfe * (c r+fc3 cos /3 fe t 



c T+k4 sin , 
Qfc * (c r+fc3 sin + c r+fc4 cos f3 k t) , 



o 

r+k 3 
~0 

y r+fc 4 
fc = 1, 



Usi(i) = ^ cosh(|w|i) + sinh(|u|i), 
v S2 {t) = v° S2 cosh(|w|i) + v° si smh.(\u\t), 

/fii(t) = ^ 2 sin(|u|t) + 5? cos(Ht), 
\v j2 (t) = ?)° 2 cos(|u|t) - v°sm(\u\t), 

J = a + 1, • • ■ , 



Cr+ki 
-k 2 j 

~ Cj— ffc3 
~ C r -\- ki , 

n 



= ^+/c 2 Sinh Sil1 + ^r+fci COsh COS M 

+Vr +kl sinh afci cos (3 k t + £° +fc2 cosh o^i sin (3 k t, 
= -Vr +kl sinh a k t sin /3 fc t + v T ° +fc2 cosh a k t cos 



^ u r+k 2 



sinh a^i cos 



Vr +kl cosh a k t sin f3 k t, 



v r +k 3 (t) = v° r+kA sinh a k t sin /3 fe i + v^. +ks cosh a fc t cos f3 k t 

-Vr +kz sinh a fc t cos /3 k t - ?)° +fc4 cosh a fc i sin f3 k t, 
v r+ k 4 (t) = -v®. +k3 sinh a k t sin /3 fc t + v® +ki cosh a fc i cos (3 k t 



+fc4 sinh afct cos /3 k t + £°+A;3 cosh a ^ sin A* 



n 
4 



Finally, we obtain the solution v(t) to the system (4.18) 
(4.19) 



sinh(\u\t) + ij° S2 cosh(|w|i) 



i (*5° 2 sinh(|u|t) + d° Sl cosh(|u|t) 



S2/ ' 



u s 1 



(4.20) 



= sin(|^|*) - S£ cos(|n|t) + 

>(*) = sin (M*) + «°i cos(|n|t) - 5° ), 



j = s + l, 



(4.21) 



A. 



u r +fci (i) = ^4^ sinh sin + cosh a k t cos 

sinh a k t cos /^i + A k cosh a^i sin /3fct — B k , 

v r+k2 (t) = —B k sinh a k t sin + A k cosh o^t cos f3 k t 

+A k sinh afct cos f3 k t — B k cosh sin f3 k t — A k 

= v^{a k v° r+k2 + (3 k Vr +kl ), B k = — 



\u\ 



u 



z{a k v Q r+kl -I5 k v r+k2 ), 



1(1 
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(4.22) 



' v r+ k 3 (t) = — Ck sinh a>kt sin fat — Dk cosh a\~t cos fat 

+Dk sinh a^t cos fat + Ck cosh a^t sin fat + Dfc, 
u r+ fc 4 (i) = Dk sinh a^i sin fat — Ck cosh a^t cos fat 

+Ck sinh a^i cos fat — Dk cosh akt sin + C&, 



Ck 



1 



1 



fc = l,... 



n r 
4 ~ 2" 



A curve £(t) has the same causal type as the initial velocity vector tr ([!])■ 

Lemma 10. T/ie projection of the curve v(t) onto the (v Sl ,v S2 ) -plane, is a brunch of the 
hyperbola with the canonical equation 



(4.23) 



v 8l + 



\u\ 



Vs 2 + 



%V»_ -(%) ' + (%) " 
til 2 



Proof. From (4.19) we calculate that VgAf) — u 2 2 (f) 
expression can be rewritten as stated in ( |4.23 ). 



^y^ ^i^-D. Th, 

□ 



Lemma 11. The projection of the curve v(t) onto the (dj 1 ,Vj 2 ) -plane, j = s + 1, . . . , r, is a 

J(i° ) 2 +(v° ) 2 

circle with the center ( — , ttt ) o/ radius — — — — — 



Proof From ( |4.20| ) we get u? (i) + 5j 2 (t) 



N 

^~ J1 ^| M |2^ sinh 2 ^ . This expression leads to 



^ 2 



+ u 



'Ja 



^ 2 _(^) 2 + (^) 2 



□ 



Lemma 12. TTie projection of the curve v(t) onto the (v r +ku Vj-^^-plane, k = 1, . . . , ? 
is a logarithmic spiral with the equation 



r 

4 2' 



(4.24) 



1 



ejf +fcl (t) +% +ka (t) = r%( (i° +fcl ) 2 + a o +fc a) 2 y ai ' 



Proof. From (4.21) we obtain 



v 2 r+kl (t) + vl +k2 {t) = W {A\ + Bl)e^. 



Substituting Ak and Bk into A k + B k we get the desired formula (4.24). 



□ 



Analogously using (4.22) can be shown the next lemma. 

Lemma 13. The projection of the curve v(t) onto the (v r+ k 3 , v r+ k 4 ) -plane, k = 1, ... 
is a logarithmic spiral with the equation 



n r 

' 4 2' 



(4.25) 



Vr+k 3 (t) + V 2 r+kA {t) 



1 



r,i° ^ 2 -HV, -2a fc t 
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4.4. Solution in a matrix form. Let us rewrite the matrix D (4.13 )-(4.15 ) in the following 
form D = PDP 1 , where the matrix P is obtained from the matrix P by rearranging lines in 
the blocks of type (4.15) and 



(4.26) 



D 



o 



V 





\u\D 2 








\u\ 2 D 3 




where D\ is a (2s x 2s)-matrix ( ^ J V D2 is a block-diagonal (2(r — s) x 2(r — s))-matrix 



consisting of (r — s) (2x 2)-blocks of type 



1 
-1 



Z?3 is a block-diagonal ((^ 



-rjx^-r 



))- 



COS C/fe 8111 ffc 

— sin 0^ cos 



where cos Ok 



matrix consisting of j — | (2 x 2)-blocks of type - 

t^p, fc = 1, . . . , ^ — |, and 1)4 = — -D3. Notice that there is no matrix D\ in (4.26) if s = 0, 
i. e. if p is even. Analogously, there are no matrices D\ and D2 in (4.26) if r = 0, and there 
are no matrices -D3 and D4 if /c = 0. 
Then 



(4.27) 



( e ^ Dit 

e \u\D 2 t 


V 









\ 



e W\ 2 D±t J 
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D l \D„ I -odd, ^ 
\ I 2s , I - even, 



D 2 , l = 4d + l, 



'2(r-s)> 



I = 4d + 2, 



D 



3,k 



D 



4,fc 



-D 2 , / = 4d + 3, 

2 (r-s), l = 4d, d = 0,1,2,... 

cos 19 k sin Wk 
- sin 19 k cos 19k 

D l k-1 - - r 



where I 2s and I^lrs) are identity matrices of dimensions (2s x 2s) and 2(r — s) x 2(r — s) 
respectively. 

Taking into account the formulae above we calculate each of the blocks in ( |4.27 ) 



\u\Dt , D 2 \u\ 2 t 2 D 3 \u\H 3 

e \u\ Dl t = hg + + — + ^LJ — + . . . 



2! 



3! 



U2 S H ^ 1 71 h...j + [\u\tDi-\ 1 \-...) 



2! 



4! y Vl 1 1 3! 

= cosh (|ii|i)/ 2 s + sinh 



i«m„f T „ i > Dl\u\ 2 t 2 D%\u\ 3 t 3 

e \u\D 2 t = j^_ s) + + 2l_J + 2l_j + 



2! 



3! 



( J 2(r- S ) 2\ + M • • ' > + y\ U \ W2 *i + ^ ■■■) 



4! 



3! 



5! 



cos (\u\t)l 2 ( r - s ) + sin (\u\t)D 2 



/ oo oo \ 

1=0 ' 1=0 

oo oo 

\ «=0 1=0 I 

e \u\H co ^ h cos (\ u \2 t sm9 k ) e\ u \ Hcosdk sin(|u| 2 t sin^ fc ) 

_ e \u\Hcose k S i n (| n |2i sin( 9 fc ) e \u\H cos e k cos (| n |2 t sin6 y 

,,,/ i cos 4/ sill . V 

- sin cos /3fci 



cos (3kt sin 
— sin /3kt cos 



fc = l,..., 



■n 



Therefore, e' u ' 2£>3 * is a block-diagonal matrix consisting of § — § (2 x 2)-blocks of type e'"' 2£>3 fc< 

and e' n ' D4i is a block-diagonal matrix consisting of | — § (2 x 2)-blocks of type e'"' 2D4 ' fct . 
The solution of the system 



(4.28) 



v = Dv, 



where v = Pv, with the initial conditions v(0) = 0, v(0) = vq, is of the form v(t) = le" 

;DK DD 1 = &D = \u\ 2 I v and D\ = D 1 ,D t 2 = -D 2 



I v )b~ 1 ^. Notice that D' 1 



l rit 
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Therefore, the solution of (4.28) reduces to 

1 



(4.29) 



u 



e Dt -I v )D%. 



In order to write the solution more explicitely, let us make the following notations: let 
v (t) denote the 2s-dimensional vector consisting of 2s first coordinates from the vector v(t), 
analogously let v 2 (t) denote the vector consisting of next 2(r — s) coordinates from the vector 
v(t), let us denote also by v 3 ' k (t) the 2-dimensional part of vector v(t) corresponding to each 
block D^k, and analogously, v A,k (t) - part of v(t) corresponding to block D^. In analogous 

way, let Vq be the 2s-dimensional vector of first 2s coordinates of the vector vq, v 2 be the 
vector of next 2(r — s) coordinates of the vector vq and let us denote by v^ k the part of vector 
vq corresponding to each block D^^, and analogously, v 
§. Then the solution Q4.29[ ) 



=4,fc 



Dd k , k = 1, 



n 



to the system (4.28) can be written in the form 

5 (*) = F72 (sinh(|n|t)/ 2s + (cosh(|u|t) - l)Di), 



part of vq corresponding to block 

=2 =3,1 =4,1 = 3i f-|^4, ? -^ (t) 



(4.30) 



sin(|n|t)/ 2(r ._ s ) + (1 - cos(|u|t))L> 2 ), 



=3,fe 



=4,fe 



where matrix 



i 

"I" 



(t) 

it) 
1 - 



=3,fc 2a fe t| 



=4,fc„-2a fc tTi 



- cos((3 k t) 
sin((3 k t) 



1 - cos(/3 fc t) 



, fc = l, 



2' 



For the horizontal part of the geodesic j(t) we make following observations: 2-dimensional 
v 1 ^), corresponding to the cell in D of type D\, describes the geodesic which lives in a sub- 
Lorentzian 3-dimensional Heisenberg group, i. e. Heisenberg group with the metric tensor 

q ^ ~\ on the horizontal bundle; 2-dimensional projections of v 2 (t), corresponding to the 

cells in D of type D2, describe the geodesies which live in a sub-Riemannian 3-dimensional 

10 



Heisenberg group with a Riemannian metric 



1 



on the horizontal bundle. There are 



no analogues to v s,k (t) or v 3 ' k (t) in Heisenberg group, but together (v 3,k (t), v^ k {t)) for each 
k they correspond to the geodesic which lives in a sub-semi-Riemannian H-type Quaternion 
group with the metric of index 2 on the horizontal subbundle and corresponds to the (4 x 4)- 
dimensional block of type 



/ cos 6k 
— sin 6k 





sin 6k 
cos 6 k 










- cos 6k 
sin 6k 



\ 


sin 6k 
cos 6 k ) 



(see [33]). 

Once v(t) is calculated one can find the horizontal part v(t) of the geodesic j(t) = (v(t),u(t)) 

by 

(4.31) v = P~ X P~ X %. 
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To find a vertical part u(t) we substitute v = P P x v into the condition (4.4) and obtain 
that 



and, integrating, we get 
(4.32) u(t) 



[P~ l p- l v, P- L p- L v] dt + u t. 



5. Example 

In this section we provide an illustrating example of how the theory above can be applied 
to concrete situation. 

Let us consider the .ff-type group with algebra V ®U , y = M 8 , U = M 7 , and with nonde- 
generate metric rj of arbitrary index p ^ 4 on V. 

Let 7(i) = (y(t),u(t)) be a geodesic on such group. It is sufficient for us to find only the 



part v(t) since v(t) and u(t) can be obtained from it by (4.31) and (4.32). 



To escape extra indexes it would be convinient to denote initial vertical velocity u by 
u = (ui, . . . ,u 7 ). Since j(u) appears from the structure of the algebra and does not depend 
on the metric on V, we will base on the fact that matrix j(u) for octonion H-type group in 
sub-Riemannian case is already calculated in [3]. 



where u 
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-u 2 




-u 3 





Ul 


Uf, 


U 7 




— U/l 


-us 
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7 . Matrices 


j a are taken from [5j 


and have the following fc 
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Let ji, . . . , jV be endomorphisms from subsection 4.2 satisfying (4.16) and such that 

j(u) = mji + . . . + u 7 j 7 . 



u\+u 2 2 +ul+u\+ul+ul+u^. 



Denote as before the matrix rjj(u) by symbol A. Denote also \u\ 2 

Let us find the characteristic polynomials and eigenvalues of A for all p = 1, . . . , 4. 
If p = 1, then 



p\(\) = A 8 + 2|u| 2 A 6 - 2|u| 6 A 2 



u 



and eigenvalues of A are ±|it|, zti|tt| , ±i|tt|, ±i|n|. Notice that all matrices rjji, 
the same eigenvalues ±1, ±i, H, ±i. 

We see that if p = 1, then s = 1, r 



, r]j 7 have 



4 and k = in (4.13)-(4.15) and, therefore, 



= (S 1 ®,^)), 5(0) = By ((430]) we get 



~v\t) = 
~v\t) 



u 



- ( sinh |ii|iT 2 + (cosh|n|t - l)£>i), 



v? 



j2 ( sin \u\Hq + (1 — cos \u\t)D2) ■ 



where D\,D2 are (6 x 6)-matrix described in (4.26), 12,16 
Up 



-- 2, then 

p* (A) = A 8 + 4n?A 6 + 2\u\ 2 (a 2 + 2{uj 



identity matrices. 



i ?))A 4 + 4|n|VA 2 + |u| 8 , 



and eigenvalues are ±(iui ± y/u 2 , + . . . + u 2 ), ±i\u\, ±i\u\. Notice that matrix r\j\ has eigen- 
values ±i, dbi, ±i, ±i and matrices r}j2, ■ ■ ■ rjjj have eigenvalues ±1, ±1, ±i, ±i. 

We see that if p = 2, then s = 0, r = 2 and k = 1, a\ = y / u|~+777+ u 2 , f3\ = u\ in 



pl3| )- pl5l ) and, therefore, 5(t) = (C 2 (t),6 3,1 (t),C 4,1 (t)) I 5(0) 



: ,2 -3,1 -4,1 



) . By ( [4301 we 
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2,'Ct 



t, 2 (t) = |^(sin|«|t/ 4 + (1 - cos \u\t)D 2 ), 



~v 3 >\t) 



=3,1 2aitT, 



1; 



4,1 -2aitT[ 



where -D2 are (4 x 4)-matrix described in (4.26), I4 is (4 x 4)-identity matrix and 

(5.1) El : 



1 / 1 — cos/3it sin fat 



If p = 3, then 
p^(A) = A 8 + 2(u? + t*l + «| - «| 



sin 1 — cos /Sit 



it^)A 6 - 2\u\ A {u\ + u\ + u\-u\ 



and eigenvalues are ±|w|, +i|u|, +(i\A*l + ^2 + u i ^ -\/ u 4 + • • • + Notice that i]ji,rjj2, rjj^ 
have eigenvalues ±1, ±i, ±i, ±i, and matrices 77J4, . . . , 77J7 have eigenvalues ±i, ±1, ±1, ±1. 
We see that if p = 3, then s = 1, r = 2 and = 1, ai = yii| + . . . + u 2 , fa = \Jv\ + it 2 , + 1X3, 

and 5(t) = (^(t),^),^)^ 4 ' 1 ^)), C(0) = (4 1 ,4 2 ^o' 1 ^o' 1 )- Then, by ((430) 



r,l 



sinh \u\tI2 + (cosh |u|t — l).Di) , 



i 2 (t) = (sin|it|tJ 2 + (1 - cos|n|t)L> 2 ), 



5 3 - 1 (t) =^' 1 e 2ait Bi, 



^(t) = ^e 



7 o 

4,1 -2aitTi 



where Bj has the same form as in (5.1) with /?i = a/u 2 + u 2 , + u 2 , D\ and -D2 are (2 x 2)- 



matrices represented in (4.26). 
If p = 4, then 



p^(A) = A 8 + 4(u 2 + u 2 + u 



+2(\u\ 4 + 2(ul + ul + u 2 3 



in 



u 4 

4/„.2 



2\ \6 



u z 7 )\ 



2 ) 2 )A 4 + 4|u| 4 (u 2 + n 2 + n 2 



Ua 



u 2 )\ 2 + \u\ 



<■!)■ 



and eigenvalues are ±(iy/u1 + u 2 , + u 2 ±^u 2 + . . . + u 2 ), ±(i-y/it 2 + u 2 , + uI + a/u 2 + • • • + 
Notice that matrices rjj± , 77/72 , have only purely imaginary eigenvalues ±i,±i,±i,±i, and 
matrices 77J4, . . . , ry'7 have only real eigenvalues ±1,±1,±1,±1. 

We see that is p = 4, then s = 0, r = 0, k = 2 and ai = a 2 = \At 4 + • • • + u 2 , fa = 

fa = \/u{ + it| + u\. Then d(t) = (6 3,1 (t), S 4,1 (t), 5 3 ' 2 (t), ?) 4 ' 2 (t)) and initial velocity is 5(0) = 
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'=3,1 z4,l =3,2 -4,2 
< v ' v > V 



). By formulae (4.30) we have 



~^\t) - 

~v^(t) - 



-3,1 2aiin 







V, 







4,2 e -2a lt j 



In fact, the author's guess is that for all sub-semi-Riemannian H-type groups satisfying j 2 - 
condition all complex eigenvalues at + iflk are the same VA; for every fixed index p. Notice that 
in case p is even there are no real eigenvalues, i. e. s = in (4.13), and in case p is odd ±|u| 
are the only couple of real eigenvalues, i. e. s = 1. 



6. Possible open problems 

As a continuation of work in this direction it would be interesting to consider the general- 
ization of matrices A to matrices of order n = p a k, where p is a prime number and k is prime 
to p, and to consider the corresponding analogue of -ff-type groups. In this case each matrix j a 
would represent a p-th root of —1 (see [IS])- The number of such matrices in a set of n-rowed 
matrices that satisfy the conditions 

3i = -I, Jaj/3 = -j/33a, (a//3), 

is not more than 2a + 1. 

Another interesting research line could be considering vector spaces U and V both with non- 
degenerate metrics (•,•)[/ and (v)v- The description of corresponding H-type groups would 
be probably simpler than in the present article, since the construction of representations of 
Clifford algebras C\(U) are less tricky in the case of nondegenerate indefinite product on U. 
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